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We report on our recent rigorous construction of complete renormalization group trajec- 
fH , tories between two fixed points for the three-dimensional phi-four model with modified 

propagator considered by Brydges, Mitter and Scoppola (BMS). These are discrete criti- 
cal trajectories which connect the ultraviolet Gaussian fixed point to the nontrivial BMS 
infrared fixed point which is an analogue of the Wilson-Fisher fixed point. The renor- 
malization group map is defined rigorously and nonpcrturbatively, without using the 
hierarchical approximation. The trajectories are constructed by a fixed point argument 
in a suitable Banach space of sequences, where one perturbs a nonlinear one-dimensional 
^ ' iteration. 

m 

, Keywords: renormalization group, crossover, heterclinic orbit, invariant manifolds, mass- 

f^*^ ■ less flow, Wilson-Fisher fixed point. 

1. Global Dynamics of Wilson's Renormalization Group 

A central mathematical problem posed by quantum field theory (QFT) is that of 
giving a rigorous meaning to functional integrals such as 



D(j). ■ . e - I ti d[^(V<j>) 2 (x) + ucl>(x) 2 +g4,( X ) 4 ]dx _ 

In this example corresponding to the so-called 4> A theory, T is the infinite- 
dimensional space of functions form M. d to K, Dcf> is the Lebesgue measure on J 7 , 
and the dots stand for some observable or functional of the field </>. A natural ap- 
proach to this problem is to try to construct such a measure on functions or rather 
distributions in iS'(R d ) by a scaling limit of lattice theories on (aL) d C R d . This 
corresponds to imposing a cut-off - on momenta in Fourier space. One can easily 
rescale to a unit lattice and see the approximants to the continuum theory as living 
on the lattice Z d instead of {aL) d . These approximants become points in the space 
of all possible unit cut-off theories. In essence, Wilson's renormalization group (RG) 
is a dynamical system acting on this space. If dv is a measure on random fields (f> 
with 4>(p) — if \p\ > 1, i.e. with sharp unit cut-off on momenta, one introduces 
a magnification ratio L > 1 and splits the field as <f> = £ + '/'low- The fluctuation 
field C corresponds to keeping Fourier modes with L^ 1 < \p\ < 1, whereas the low 
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momentum or background field 4>\ ow has Fourier modes with \p\ < L . One can 
then integrate over £, and obtain the marginal probability distribution on </>i ow com- 
ing from the original probability distribution dv. Finally one rescales back to unit 
lattice by letting ip(x) — L^<f>\ OVJ (Lx), and considering the resulting probability 
distribution dv' on ip. The RG map is the transformation dv — > dv' . Important 
features of this dynamical system are the notions of fixed points, eigenvalues of lin- 
earized RG around them, local stable and unstable manifolds etc. They play a key 
role in controlling scaling limits and constructing continuum functional integrals 
by lattice approximations. The features just listed which have also been the most 
studied ones are local, i.e., pertaining to the vicinity of single fixed points. Much 
less is known about global features of the RG dynamical system, e.g., heteroclinic 
trajectories between fixed points. This is the main object of this article which briefly 
reports on recent progress we made in this area. 1 

2. Rigorous Results (Selection) 

There is a long history of rigorous results on the RG dynamical system. We will 
only mention a selection of such results which gives the best idea of the context into 
which the results of our article 1 fit. 

The RG exponents for the Gaussian fixed point, the existence of a nontrivial 
infrared fixed point in '4 — e' dimensions, as well as the construction of its local sta- 
ble manifold, were obtained by Bleher and Sinai, 2 ' 3 in the hierchical model (HM) 
approximation. Similar constructions were also provided by Collet and Eckmann, 4,5 
Gawedzki and Kupiainen. 6,7 Furthermore, the two-dimensional local unstable man- 
ifold of the Gaussian fixed point for the r/>| HM was constructed. 8 The hierarchical 
</>! infrared fixed point (i.e. for e = 1) was also constructed. 9 The appearance of 
new fixed points at dimensions d = 2 + n — 3,4, .. . in the local potential 
approximation was established by Felder. 10 

As per the Euclidean model, i.e., the full model without hierarchical or local 
potential aproximations, the existence of a nontrivial infrared fixed point as well 
as the control of its local stable manifold was provided by Brydges, Dimock and 
Hurd. 11 Later, a similar result for a nicer '(4 — e)-dimensional' model was derived 
by Brydges, Mitter and Scoppola. 12 

The previous results concern local features of the RG. Rigorous results of a 
global nature are more scarce. One can mention the proof of uniqueness of the 
infrared fixed point in the local potential approximation in dimension 3 < d < 4 
by Lima. 13 Closest to our result is the construction of the massless Gross-Neveu 
model in '2 + e' dimension by Gaw§dski and Kupiainen. 14 Their result amounts to 
controlling a heteroclinic trajectory joining a nontrivial ultraviolet fixed point to 
a trivial infrared fixed point, for a Fcrmionic QFT. In our article 1 which concerns 
the Bosonic BMS model, 12 we obtained the construction of discrete heteroclinic 
trajectories joining a Gaussian ultraviolet fixed point to a nontrivial infrared fixed 
point. One should finally also mention related work on massless nontrivial QFT 
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models such as (f> 4 in two and three dimensions 15-17 or the Thirring model. 18 
3. The BMS Model 

The BMS model for a scalar field <p ■ IR 3 — > K corresponds to the functional integral 

potential V(4>) 

Z = j D<j> e -3<^(- A )^ 1 *>I,2(IR3) -'/^(^ 4 WW(^)0 . 
V v ' 

Gaussian measure 

The free propagator is (— A)~t~ (x, y) ~ j^~pi?] where [0] = ^j 1 denotes the 
canonical scaling dimension of the field. One also has a multiscale decomposition 
of this propagator, in the form / °° ^l~ 2 m u (^j 2 ) with u a smooth, compactly 
supported, rotationally symmetric function which is nonnegative in both x and p. 
Instead of a lattice regularization, it is more convenient to use continuous cut-offs, 
by taking C(x — y) = ^l~ 2 W u (^j^) as a free propagator. The high versus 
low momentum split corresponds to the decomposition C(x — y) = T(x — y) + 
C L -i(x — y) with C L -i(x — y) = L^ 2 ^C(L^ 1 (x — y)). The fluctuation covariance 
therefore is T(x — y) = ^l~ 2 W u (^j^)- One can rewrite accordingly the original 
Gaussian measure as a convolution d^c = d^v *d\ic L _ x - The functional integral, 
with integrand Z(4>), becomes 

Z = J d» c (<P) Z(^) = J dpc^WWriQ Zty + Q = J d(, c (<P) {nz){4>) 

where {nZ){(j>) = J dfj, T (C) Z((f) L -i + () and 4> L -i(x) = L~^cj)(L~ 1 x). Now the 
RG map can be seen as acting on integrands: Z — ► 1ZZ. 



4. Good Infinite Volume Coordinates 

An important ingredient for the rigorous implementation of the RG is a good set of 
coordinates for the dynamical variable Z, which behave well in the infinite volume 
limit. Such a formalism was introduced by Brydges and Yau. 19 One uses the decom- 
position of R 3 into unit cubes or cells defined by a lattice Z 3 . If A is the large finite 
box where the model is defined, one writes the following polymer representation for 
the integrand: 

oo ^ 

Z ( A ^) = E^ E ex p 

n— X-i,...,x n 

disjoint in A 

xK(X 1 ,4>\ Xl )---K(X n ,4>\ Xn ) . (1) 

The Xi arc polymers, i.e., connected finite unions of closed unit cells. The functionals 
AT (A, (j)) are local in the sense that they only depend on the restriction of the field to 
the sets X. The coordinates for Z now are triples (g, fx, A) where g, ^ are the running 
couplings for the 4 and mass terms, and A = (A(A, -))x polymer is a collection of 



/ dx{g : (f> 4 (x) : c +n : 4> 2 (x) : c } 

JA\(UXi) 
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local functionals. The need to explicitly extract second order perturbation theory, 
in order to analyse the RG dynamics, imposes an additional rewriting of the form 

K(X,<f>) = g 2 [explicit complicated forrnula]e _v(x ' 0) + R(X,<f>) . 

So the dynamical variable becomes the collection of functionals R which can be 
thought of as an 0(g 3 ) remainder. As a result, the RG for the BMS model can be 
written as an explicit (albeit complicated) map (g,fi,R) — ► (</, //, R'). It has the 
form 

g' = Ug - L 2 *a(L, e)g 2 + £ g (g, fi, R), 
H 1 = L^fi + £^{g, fx, R) , 
R' = £^(R) + Z R (g,»,R) . 

The quantity a(L,e) is ~ log£. The £'s are small nonlinear remainder terms. The 
C 9 '^ stands for a (g, /z)-dependent linear operator in R. Its crucial property is that 
it is contractive, which follows from the so-called "extraction step" in the definition 
of the RG map. This is the analogue in the present setting of the traditional BPHZ 
subtractions in renormalization theory. 20 It removes the dangerous local part of 
the two and four point dependence in the low momentum field cj>^-i, from the 
functionals R and transfers them to the exponentiated potential. This relies on the 
nonuniqueness of the polymer representation of Eq. [T]for a given integrand Z. 

Ignoring the £ remainders makes for a simplified RG map with easy dynamics. 
In addition to the trivial fixed point (0,0,0), it has a fixed point at (g*,0,0) with 
9* = T^a ^ e - It also has heteroclinic trajectories (g n , 0, 0)„ e z parametrized by the 
single number go in the interval (0, g*), and constructed using the one-dimensional 
iteration by the function f(x) = L e x — L 2e ax 2 . The BMS fixed point was obtained 
as a perturbation of the approximate fixed point (<7*,0,0). Similarly, the true RG 
connecting orbits are constructed by perturbing the above approximate trajectories. 
Our main result 1 is as follows. 

Theorem 4.1. If the bifurcation parameter e > is small enough, for any ujq £ 
]0, g[, there exists a (locally unique) complete trajectory (g n , fJ, n , Rn)n& 

for the RG 

map such that lim (g n ,fj, n ,R n ) — (0,0,0) the Gaussian ultraviolet fixed point, 

n — > — oo 

and lim (g n ,/i n ,R n ) = (g* , /i* , i?* ) the BMS nontrivial infrared fixed point, and 

n — >+oo 

determined by the 'initial condition' at unit scale go = ujo9*- 
5. Idea of The Proof 

One writes (g n , fi n , Rn) — (g n + Sg n , \x n , R n ), so the trajectory equations are 
rephrased in terms of the deviation variables (Sg n , \x n , R n ) relating the true tra- 
jectory to the approximate one. One has to enforce three boundary conditions: 
1) jji n does not blow up when n — * +oo (infrared), 2) R n does not blow up when 
n — > — oo (ultraviolet), 3) Sgo = (unit scale). One iterates the flow equations either 



October 6, 2009 17:41 



WSPC - Proceedings Trim Size: 9.75in x 6.5in ICMP2009 



5 



backwards or forwards towards the appropriate boundary condition. For instance, 
the fj, n term is reexpressed by forward iteration as: 



The R is iterated backwards. Finally the g n is iterated backwards if n > 0, and 
forwards if n < 0. The four 'integral' equations, including Eq. [5J obtained in this 
way are interpreted as a fixed point equation in a big Banach space of two-sided 
sequences. The sought for trajectories are thus obtained using Banach's contraction 
mapping theorem. 

6. Functional Analysis, Norms 

o 

Let A denote a (closed) unit cell. There is a Sobolev imbedding FT 4,2 (A) c — > C 2 (A). 

o 

The fields <j> live in Hilbert spaces Fld(X) = acX W 4,2 (A) where extra C 2 gluing 
conditions are imposed, using the Sobolev imbedding. The norm on a field is given by 
H^lliicUXI = ICacxzCm<4 ll'^'MII 2 o ■ The fluctuation measure d//r is realized 

in the Hilbert spaces Fld(AT). One also needs a second field norm ||0||c 2 (jf) = 

sup max \d" cj> (x) I . While the r/>'s are real- valued, the functionals can be complex 

xex kl<2 

valued and are compared using norms of the form 



where D n K denotes the n-th order differential. One of the main technical difficulties, 
or fibered norm problem, is that these norms depend on the dynamical variable g. 
The way around it is to use the approximate solution g n to calibrate them. 1 

7. Perspectives 

Many open problems in the continuation of this work remain. One should study more 
refined dynamical systems features of the BMS model such as the construction of 
the full hctcroclinic invariant curve. One should study its regularity properties and 
smoothness at g = (asked by K. Gawe,dzki). There is the problem of constructing 
the correlation functions for the primary field r/> and also for composite fields and 
study possible anomalous dimensions (there are preliminary results on this by P. 
K. Mitter). A good testing ground for these more ambitious goals is the case of a 
hierarchical version of the BMS model (ongoing work by the author together with 
G. Guadagni and Ph. D. student A. Chandra). This could lead to the construction 
of the massless Euclidean BMS model over the p-adics. If the previous hurdles are 






X 
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cleared, one should also study the analytic continuation to Minkowski space (in the 
Archimedean case). Another interesting project would be to provide an alternate 
construction of the critical BMS theories using recent phase space expansion tech- 
niques. 21 Finally the problem of controlling complete RG trajectories (from the UV 
to the IR end of the scale spectrum), here addressing the BMS model, is a worthy 
and certainly difficult task for many more QFT models: Gross-Neveu in 2d, 4>^ at 
large N , noncommutative vulcanized r/>|, the 2c? a-model, Yang-Mills in 4c?. The 
inspiring talk by S. Weinberg on QFT and the asymptotic safety scenario, at this 
congress, also leads one to wonder if one could devise a rigorous RG framework in 
order to study a nontrivial ultraviolet fixed point in quantum gravity. . . 
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